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Abstract 



The asymptotic collinear factorisation theorem, which holds for diffractive deep-inelastic 
scattering, has important modifications in the sub-asymptotic HERA regime. We use per- 



turbative QCD to quantify these modifications. The diffractive parton distributions are 
shown to satisfy an inhomogeneous evolution equation. We emphasise that it is neces- 
O sary to include both the gluonic and sea-quark t-channel components of the perturbative 

■ Pomeron. The corresponding Pomeron-to-parton splitting functions are derived in the 

Appendix. 

> 

X 1 Introduction 

a: 

A notable feature of deep-inelastic scattering is the existence of diffractive events, 7*p — > X +p, 
in which the slightly deflected proton and the cluster X of outgoing hadrons are well-separated 
in rapidity. 1 At high energies, the large rapidity gap is believed to be associated with Pomeron, 
or vacuum quantum number, exchange. Some secondary Reggeons also have vacuum quantum 
numbers, but these contributions are exponentially suppressed as a function of the gap size, 
and are negligible at small xjp. The diffractive events make up an appreciable fraction of 
all (inclusive) deep-inelastic events, 7*p — > X. We will refer to the diffractive and inclusive 
processes as DDIS and DIS respectively. 

The recent improvement in the precision of the DDIS data [1, 2] allow improved analyses 
to be performed and more reliable diffractive parton densities to be extracted. Moreover, 
measurements of diffractive charm production in DIS are available [3], which provide added 
constraints, particularly on the diffractive gluon density. 



1 The process and the corresponding variables are shown in Fig.^ 
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Figure 1: Diagram showing the kinematic variables which describe the DDIS process 7*p — > 
X + p: is the fraction of the proton's momentum transferred through the rapidity gap, 
(3 = x B /xip is the fraction of this momentum carried by the struck quark, x B is the Bjorken x 
variable, q 2 = —Q 2 is the photon virtuality, and t is the squared 4-momentum transfer. 



It is common to perform analyses of DDIS data based on two levels of factorisation. First 
the diffractive structure function F2' may be written as the convolution of the usual coefficient 
functions as in DIS with diffractive parton distribution functions (DPDFs) a D [4]: 

if (3) = £cV®a D , (1) 

with factorisation scale /x F , where aP = (3q B or (3g B satisfy DGLAP evolution in \i F . The 
collinear factorisation theorem (pQ) applies when /j, F is made large, therefore it is correct up 
to power-suppressed corrections. In a second stage, Regge factorisation is usually assumed [5] 
(see, for example, the preliminary HI analysis [2]), such that the diffractive parton densities aP 
are written as a product of the Pomeron flux factor fip{xip, t) and the Pomeron parton densities 
cF = [3q IP or ^cf '. Taking \i F = Q, the t-integrated form is 

a D (xjpJ,Q 2 ) = faM^iPiQ 2 ), (2) 
where the Pomeron flux factor is taken from Regge phenomenology, 

/^min gSjp t 

fipM = / dt 2ajp(t) _ i ; (3) 

with ajp(t) = ajp(0) + a'jpt. For simplicity of presentation we omit the contribution of secondary 
Reggeons to the right-hand side of (J2J). Strictly speaking, the parameters ajp(O), a'jp, and Bjp 
must be taken from fits to soft hadron data. However, then one fails to describe the xjp 
dependence of the diffractive data. Therefore, as a rule (see [2]), a'jp and Bp are fixed 

by the analyses of soft hadron data, while ajp(O), and the parameters describing the input 
Pomeron parton distributions at some scale Hq, are determined from a fit to the DDIS data. In 
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these analyses, the Pomeron is treated as an effective pole in the complex angular momentum 
plane, and regarded as a hadron-like object of more or less fixed size. This Regge factorisation 
takes place in the non-perturbative region at some low scale /i, with /i < /i , of the order of the 
(inverse) size of the hadron. However, in such fits, the value of ajp(0) extracted from DDIS data 
(for example, 1.17 in Ref. [2]) lies significantly above the value of 1.08 obtained from soft hadron 
data [6]. This indicates that there is a contribution coming from the small-size component of 
the Pomeron, that is, coming from the perturbative QCD region where the vacuum singularity 
has a larger intercept. 

This Regge factorisation approach is a simplified phenomenological model. Here we shall 
not assume Regge factorisation for the whole a D , but instead study the impact of applying 
perturbative QCD to the analysis of DDIS data. The procedure we discuss in this paper 
formed the basis of our earlier DDIS analyses [7,8]. 

The content of this paper is as follows. In Section |21 we recall the collinear factorisation 
which underlies the description of inclusive DIS. This paves the way for the discussion in Section 
El of the evolution equations for the diffractive parton distributions aP due to the perturbative 
Pomeron. These equations contain an appreciable inhomogeneous term 2 , analogous to the 
inhomogeneous term in the evolution equations for the parton distributions of the photon. 
Since the evolution equations for the diffractive densities (X cLFG cL little subtle, it may be helpful 
at this stage to look ahead to the discussion in Section Although collinear factorisation 
holds in DDIS in the asymptotic limit, at the relevant HERA energies there exist important 
modifications. These are discussed in SectionHJ Then in Section|S]we explain why it is necessary 
to include a Pomeron made of a sea-quark-antiquark pair in addition to the Pomeron made of 
two t-channel gluons. Section El shows some relevant properties of a recent analysis [7] using this 
approach, and Section [7| summarises how universal diffractive parton densities are extracted, 
after allowing for the modifications to DDIS factorisation. The Appendix contains a derivation 
of all of the Pomeron-to-parton leading order (LO) splitting functions that are necessary to 
analyse DDIS data. 

2 Collinear factorisation 

The problem in the analysis of both DIS and DDIS data is that we can only use perturbative 
QCD (pQCD) at small distances, that is, at large Q 2 . Within pQCD we can study the evolution 
of parton distributions, but the initial distributions at some relatively low scale fi are of non- 
perturbative origin and, at present, have to be determined by fitting to the data. A factorisation 
theorem underlies the analysis. It enables the amplitude to be factored into two parts, one 
purely in the pQCD domain, and the other parameterised by a phenomenological ansatz. In 
terms of Feynman diagrams, the factorisation is based on the resummation of the series of the 
most important higher-order corrections where the small coupling as is enhanced by a large 
\n(Q 2 /fi 2 ). That is, we can divide such diagrams, at a 'logarithmic loop or cell', into a part 
depending only on large scales from a part containing the low scale. Let us be more explicit. 

2 The existence of the inhomogeneous term has been known for some time [9, 10], but it is ignored in most 
phenomenological analyses of DDIS data. 
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First, recall how collinear factorisation occurs in DIS. In a physical gauge (such as the 
axial gauge for the gluon field with A^q'^ 1 = 0, where q' is the light-cone vector in the photon 
direction) the leading log contributions come from ladder-type graphs 3 ; see Fig. [2] [11]. In such 
graphs any box may be considered as the logarithmic loop which provides the factorisation at 

2 

some scale \i = fj, F . Indeed, the integral over each virtuality qi takes the form J qi+1 dqj/qf, and 
in order to generate a large log to compensate the small as, we need to be in the strongly- 
ordered region 4 with qf <C qf +1 . Because of this strong ordering, the parton may be considered 
on-mass-shell for the upper part of the diagram, which can be regarded as the matrix element 
of the hard subprocess where all the other virtualities are much larger than g«. This upper 
part represents the 7*^ interaction, and contains the coefficient function and the upper part 
of the ladder. At LO we keep the LO coefficient function and the LO splitting functions, so 
each upper loop is logarithmic and could provide an alternative factorisation scale fi F . The 
distribution of parton qi depends on the parameterised input distribution at some fixed scale 
/j,o determined by fitting to the data in the pQCD domain. Starting with the phenomenological 
distribution at fi the PDF is evolved up to the factorisation region \i F . The physical result 
does not depend on the choice of the factorisation scale \i F (up to 0(a s ) corrections when 
working at LO accuracy). If we change fi F we move from one cell to another. The change in 
the distribution in going to the new \i F is compensated by the change of the hard interaction 
in the upper part of the diagram. 

Corrections to this LO picture occur if two gluons may have comparable virtualities or if 
an extra gluon were to cross one or more of the existing gluon rungs (as shown in Fig. |5J). In 
these cases the large log integrations are absent. In such contributions the extra powers of as 
are unaccompanied by logarithmic enhancement, that is, they are higher-order in as (NLO, 
NNLO, . . . ) contributions to the splitting functions. If the additional gluon couples to the 
upper blob or the uppermost rung, then it would involve higher-order contributions to the 
coefficient function. 

3 Evolution of the diffractive parton densities 

Let us now turn to DDIS. Here we have to include a rapidity gap between the ladder and the 
proton target, and to identify the appropriate factorisation scale for the DPDF evolution. 

3.1 Starting scales 

Consider the LO Feynman graph where Pomeron exchange is described by a gluon ladder; see 
the left-hand side of Fig. El The kinematics of the process are fixed by the momentum k of the 
first emitted parton at the end of the rapidity gap. Indeed, the virtuality of the first t-channel 
parton in the upper ladder fixes the scale /1 2 = kf//3, where j3 is the (light-cone) momentum 

3 Bcsidcs the ladder graphs, we have to include virtual loop corrections, which may be included in the usual 
way by the plus prescription. 

4 (LO) DGLAP evolution effectively sums up these leading ln(Q 2 //i 2 ) terms. In the case of BFKL we have 
an analogous factorisation. However, instead of / dqf/qf, we have J dzi/zi to provide the strong ordering in 
emission angles of the gluons. In this case the ^-channel partons are two reggeised gluons. 
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Figure 2: Ladder-type diagram for the DIS cross section. The gluon line crossing two gluon 
rungs means that this is a NNLO contribution with two a$ factors unaccompanied by ln(Q 2 //^ 2 ) 
enhancements. 



fraction of the Pomeron carried by parton k and k t is its transverse momentum. The derivation 
of n 2 = k 2 /f3 is given in the Appendix; see (jA.44|) . We assume that \t\ <C k 2 . This scale /i plays 
the role of the lowest factorisation scale for the usual DGLAP evolution in the upper part of 
the diagram. Simultaneously, it is the upper scale for the lower part of the diagram. Indeed, 
the integral over the transverse momentum l t of the ^-channel gluon has the logarithmic form 
for I 2 <C fi 2 , whereas for I 2 3> /i 2 it converges as dl 2 /lf and its contribution may be neglected as 
usual. After accounting for the evolution in the lower part of the diagram, and integrating over 
l t , the probability amplitude to find the appropriate t-channel gluons is given, at LO, by the 
conventional 5 integrated gluon distribution of the proton, xjpg{xjp, fi 2 ), which is known from 
global analyses of DIS and related hard-scattering data. 

The evolution of aP is a little subtle. Before we integrate over the momentum k t — that is, 
selecting events with a fixed transverse momentum of the lowest parton — we can indeed see 
that the scale \i is the lowest possible factorisation scale for the diffractive parton densities aP . 
However, for inclusive DDIS, we must integrate over k t . The integral over k t translates into an 
integral over /i, and may be written, up to a normalisation factor, in the form 6 



d/ ' 2 1 " s 2 - a"'. (4) 



l> 



2 [I' Xjp 



A 4 



The term [...] 2 /xip plays the role of the Pomeron flux, which occurs in the conventional 
analyses; see (J2J) and (JHJ)- At first sight, the integral is concentrated in the infrared region of 
low fi. However, for DDIS we consider very small xjp. In this domain the gluon has a large 
anomalous dimension. Asymptotically, BFKL predicts xjpg(xjp, /i 2 ) ~ (/i 2 ) 5 for fixed as [13]. 
In this case the integral (JH) takes a logarithmic form, and we cannot neglect the large scale fi 
contribution. This complicates the evolution by adding a non-negligible inhomogeneous term 
to the DGLAP equation for the diffractive parton densities aP . 



3.2 Inhomogeneous contribution to the evolution equation for a 

Before we present the evolution equation for DDIS, it is informative to recall the evolution 
equation for the parton distributions of the photon, which also contains an inhomogeneous 
term. It is of the form [14] 

where, for simplicity, we show only the quark non-singlet evolution. The right-hand side of the 
evolution equation now contains an inhomogeneous term, which acts as an extra source of quarks 
produced by the splitting of a point-like photon. The point-like photon has a distribution of 
the form 

l(y,Q 2 ) = 6(l-y). (6) 

5 To be precise, the skewed gluon distribution is required, which can be written as a constant factor l|A.23|) 
multiplied by the conventional gluon distribution [12]. 

6 The derivation of the form of 10} is given in Section TA. II of the Appendix. 
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Figure 3: A ladder-type diagram showing a contribution to the diffractive parton densities 
a D (xip, /?, fj,p) in the perturbative region, fi > /i ~ 1 GeV. DGLAP evolution for the Pomeron 
densities (aF = (3(f ', fig^) is performed from fi to fiF for each component \i of the perturbative 
Pomeron in the perturbative interval \x§ < \i < \ip and then the sum is taken. The virtuality 
/i 2 of the first t-channel parton in the upper ladder is fixed by the momentum k of the emitted 
parton at the edge of the rapidity gap, /i 2 = k\ / /?; see (IA.44J) . The dashed lines in the upper 
ladder may be either gluons or quarks. The two gluon ladders shown in the diagram represent 
Pomerons. In general, these ladders can also contain quarks as well as gluons. We label the 
contributions where the two uppermost t-channel partons (labelled l± on the left-hand side) are 
gluons, in both lower ladders, by IP = G. There are also contributions in which both ladders 
have a sea-quark-antiquark pair as the two uppermost t-channel partons. We denote these 
contributions by IP = S. The flux factor in this case, fp=s, is given by ((7J) with g replaced by 
the sea-quark density of the proton, S. There are also interference terms in which one ladder 
is IP = G and the other ladder is IP = S. We denote this contribution by IP = GS. The 
factorisation scale \j? F is usually taken to be Q 2 . 
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Note that, since the photon is point-like, the extra source is independent of Q 2 . 

We also have an extra source of contributions in the evolution of the diffractive parton 
densities, a D , for ix in the perturbative region. Thus the evolution equation for aP contains an 
inhomogeneous term. This term depends on the scale \i. From (J1J, we see that the dependence 
is described by the factor 



fp(xip] /i 2 ) = 



xjp I n 



■X]pg(x]p,iJ, J 



(7) 



which specifies the strength of the inhomogeneous source of partons coming from the component 
of the Pomeron of size (Our choice of the normalisation factor N is specified in the 

Appendix; see ()A.30|) .) To be more precise, we mean partons coming from the component of 
the Pomeron wave function corresponding to the integration of the lower parts of the diagram 
shown in Fig. |3] over l t up to scale ix. An explicit pQCD calculation [15] shows that the 
corresponding splitting functions of this perturbative Pomeron into quarks and gluons are of 
the form 

P qlP (P) oc (3 3 (1-P), (8) 
P gIP (f3) oc (l + 2/3) 2 (l-/3) 2 . (9) 

The derivation and normalisation 7 of the splitting functions (and flux factor fjp) are given in 
the Appendix; see (jA.31|) and (jA.58|) . Differentiating (jlj) with respect to InQ 2 , we see that the 
evolution equations for the diffractive parton densities (a D = /3g D , /3g B ) are 

doP [ Q2 d/j 2 dcF 

fjp(xjp; /x 2 ) + f^x^^aFi^Q 2 ]^ 



d\nQ 2 J# ^ J ^^ J d\nQ 2 ' 



/IT < * 



^ V Paa>®a' D + Mxjp-Q^Pojp^). (10) 

/.IT < J 



2tt 



a=g,g 



Here, aF^, Q 2 ; fi 2 ) are the Pomeron parton densities DGLAP-evolved from a starting scale /i 2 
up to Q 2 , from input distributions oF(f3, /z 2 ; /x 2 ) = PaipiP). 

If we were to assume that fjp were independent of /1 2 , as is the case for BFKL asymptotics, 
then ()10|) would be an inhomogeneous DGLAP equation exactly analogous to that for the 
photon. The only essential difference would be the form of the j3 dependence of the splitting 
functions Paip- At first sight this appears strange. The Pomeron, unlike the photon, is an 
extended object, and we might have anticipated a form factor dependence fjp ~ l/(fi 2 Rjp). For 
such an extreme behaviour, the inhomogeneous term would be a power correction and could 
perhaps be neglected in the evolution equation. However, as xjp — > 0, the proton gluon density 
g(xip, fi 2 ) grows as (fi 2 ) ' 5 , and compensates the form factor power-like suppression. The HERA 
domain is an intermediate region: the anomalous dimension 7 is not small, but is less than 0.5, 

7 Of course, only the product of the flux factor and the splitting functions have a precisely defined normali- 
sation. Our choice for the normalisation of the separate factors is given in the Appendix. 
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Figure 4: The /i 2 dependence of the flux factor fjp, given by (jJJ) with N — 1, for three different 
values of xjp and using the MRST2001 NLO gluon distribution of the proton [16]. 

with g ~ (/W 2 ) 7 - Thus, although the integral in (@J) is convergent at large /i 2 , we still cannot 
neglect the inhomogeneous term in (|TU|) . 

The behaviour of the inhomogeneous term can be seen from Fig. which shows the /i 2 
dependence of the flux fjp of (JTj) with iV = 1 for three typical HERA values of xjp. Of course, 
Jjp is not constant and dies out with increasing ff 2 , but it is still appreciable up to rather large 
/i 2 . Therefore, we cannot neglect the inhomogeneous contribution when we describe HERA 
diffractive data with Q 2 ~ 10-50 GeV 2 . Typically the perturbative contribution, from scales 
H > hq ~ 1 GeV, is found [7] to be responsible for about half or more of the DDIS HERA 
data, depending on the value of (3. Moreover, it is seen from Fig. HJthat, even for very large 
Q 2 , DDIS factorisation must be modified when working with input scales /i 2 , less than about 
10-50 GeV 2 . Also, note that the convergence of the /i 2 integral decreases with decreasing xjp. 

In practice, rather than to solve the inhomogeneous equation (J 10)) directly, it is more con- 
venient to use a separate standard homogeneous DGLAP equation for each component fi of 
the Pomeron. Then we can sum up (that is, integrate over /z) the different /i 2 contributions, in 
which the standard DGLAP evolution of each component starts from its own scale /x, provided 
that /i is in the pQCD domain (fi > /i ), and continues up to the (collinear) factorisation scale 
fip- The contribution coming from fi < fi must be treated non-perturbatively. This component 
of the diffractive densities aP is included in the starting distributions at /io whose parameters 
are obtained by fitting to data. 



9 



4 Factorisation in DDIS and sub-asymptotic effects 



From the formal viewpoint, for any fixed xjp, the inhomogeneous term in the evolution equation 
f)10jl for aP dies out as Q 2 — > oo. Thus if we evolve from a large enough starting scale, /is, we 
are entirely in the perturbative regime and the diffractive parton densities, a D , are described 
by the usual (homogeneous) DGLAP equations and satisfy the collinear factorisation theorem, 
just like the parton densities for DIS whose factorisation was described in Section |21 The 
factorisation theorem for DDIS was proved in Ref. [4]. 

However, this is not true in the HERA domain, and in this section we emphasise those 
places where a more precise treatment of diffractive parton densities is necessary. Neverthe- 
less, universal diffractive parton distributions can still be obtained from analysing DDIS data, 
provided that care is taken in the analysis. 

4.1 Inclusion of the inhomogeneous term and the effect on g B 

For a low starting scale /is, but still satisfying /is 3> Aq CD , we can no longer omit the inhomo- 
geneous term in the evolution ()10|) of the diffractive parton densities, aP . One consequence is 
that we will obtain a smaller diffractive gluon density g D . Indeed, it is known from the global 
analyses of DIS data that the gluon is mainly driven by di^/dlnQ 2 ~ otsg- However, for 
DDIS, the perturbative Pomeron contribution to is [7] 

r Q2 du 2 

/ ^Mx^/i^Ff^Q 2 ;/! 2 ) (11) 

^ dF "d^ 2) + M** Q 2 ) F *^ G 2 ; Q 2 )- ( 12 ) 

Here, _F,f (/3, Q 2 ; /i 2 ) is the Pomeron structure function, evaluated from the DGLAP-evolved 
Pomeron parton densities a iP (/3, Q 2 \ /i 2 ). The first term of (fH?J) behaves roughly as asg D - 
Therefore, part of the derivative <9-F 2 D( ^/<91n(5 2 comes from the upper limit Q 2 of the integral 
over the Pomeron scale /z, and so this results in a smaller diffractive gluon density than if the 
second term of (|T2j) was neglected. 

4.2 Heavy quark contributions 

To be specific we consider the contributions of the charm quark, c. For moderate Q 2 it is 
convenient to use the fixed flavour number scheme, where the charm contribution arises from 
photon-gluon fusion 7*5^ — > cc. However, in this case we will miss the diagram where the 
Pomeron directly produces the c quark, that is, when c is the lowest parton k in the upper 
ladder of Fig. El In this case there is no evolution and no factorisation. This contribution 
should be added separately. 

At high Q 2 the difference between light and heavy quarks disappears. Then we should use a 
variable flavour number scheme, where a diffractive charm density is introduced at scales above 
the charm threshold. 



? D(3) 
2, pert. 



{X!P,P,Q Z 



UJr 2,pert. 

«91nQ 2 
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Here is a good place to emphasise the difference between the Pomeron densities a and the 
t-channel content of the Pomeron. Already in the charm example above, we see that although 
the Pomeron is built up only of gluons, it produces not only gluons, but also light and heavy 
quarks. 

In analogy to heavy quark production, the direct coupling of the perturbative Pomeron has 
to be included in the description of high E T dijets in DDIS. Again, this cannot be described 
purely as a convolution of a D with a 'hard' matrix element. This direct contribution should be 
added separately. 

4.3 Twist-four contributions in DDIS 

From (jSJ) and (JUJ), we see that the leading-twist splitting functions vanish as (5 — > 1. For large f3 
and the Q 2 values typical at HERA, we cannot neglect the twist-four contribution, which 
goes to a constant value as /3 — > 1. This contribution was calculated in Ref. [15], and turns out 
to be numerically appreciable. Moreover, it dies out with increasing Q 2 rather slowly since the 
extra 1/Q 2 twist-four suppression is partly compensated by the sizeable anomalous dimension 
7 of the gluon; recall [g(xjp,Q 2 )] 2 ~ (Q 2 ) 2,1 ■ This twist-four contribution goes beyond the 
leading twist-two approximation. It is described by its own twist-four evolution and coefficient 
functions. 

5 The sea-quark component of the Pomeron 

So far we have taken the Pomeron to be a parton ladder where the two uppermost t-channel 
partons are gluons. However, at small scales /i 2 , the gluon densities have a valence-like structure. 
They decrease with decreasing x already from x ~ 0.01 for /i 2 ~ 2 GeV 2 ; see Fig. |S| On the 
other hand, the sea-quark density S = 2(u + d + s) increases as xS ~ x -0 ' 2 with decreasing 
x. As a consequence we have to include another contribution to the Pomeron in which the 
two uppermost t-channel partons in the lower ladders in Fig. |H] are a sea-quark-antiquark pair. 
Thus we must introduce a sea-quark Pomeron flux given by (J7J) with xjpg replaced by 

xipS. To avoid confusion we denote the flux in (J7J) by fip=G- In general, the two lower ladders 
in Fig. 01 shown as gluon ladders, contain both quarks and gluons. For a gluonic Pomeron 
(IP = G) or a sea-quark Pomeron (IP = S) the two uppermost partons in the ladders are 
gluons or sea quarks respectively. In addition, we must include the interference contribution 
coming from one IP = G ladder and one IP = S ladder in Fig|21 We denote this contribution 
by IP = GS. 

We present the calculation of all the LO Pomeron-to-parton splitting functions in the Ap- 
pendix. The calculation of the splitting functions for the gluonic Pomeron (see P q jp=G an d 
P g jP=G of (jSJ) and respectively) already exist 8 in the literature [15], but those for the sea- 
8 Here we clarify certain factors of 2. 



11 



~ I 1 — I — I I I 



T 



~i — i — i i 1 1 1 



~i — i — i — i i i 1 1 1 



~ i — i — i — i i i 



Q 2 = 2 GeV 2 



MRST2004 NLO 



CTEQ6.1M 




Figure 5: The behaviour of the gluon and sea-quark distributions at Q 2 = 2 GeV 2 found in the 
CTEQ6.1M [17] and MRST2004 NLO [18] global analyses. The valence-like behaviour of the 
gluon is evident. 



quark Pomeron are derived for the first time. They have the forms 

P q jp=s(P) oc (3(1- (3), 
P g p=s{P) oc (1-/3) 2 ; 



(13) 
(14) 



see flA.71|) and (jA.84J) of the Appendix. In addition, there are interference contributions between 
the gluonic Pomeron and the sea-quark Pomeron, which we label by the notation IP = GS. 
These splitting functions have the forms 



gJP=GS 

see flA~89l) and (lA~96l) of the Appendix. 



P q p=Gs(j3) oc (3 2 (l-f3), 
P n jP=rAP) oc (l- / 3) 2 (l + 2/?); 



(15) 
(16) 



6 pQCD analysis of DDIS data 

The pQCD approach described above has been used to analyse the new HERA DDIS data [7]. 
A good description of both the ZEUS and HI data was obtained. The plots in Fig. |H1 show the 
separate non-perturbative and perturbative (IP = G,S, and the GS interference) contributions 
to F2^ 3 \ for two sets of values of xjp, /3, and Q 2 . The numbers in brackets in Fig. |U1 show 

9 These forms were used in the analyses of Ref. [7, 8] . 
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the total contribution of the different components starting from fi = 1 GeV; and correspond 
in the perturbative cases to taking the integral over /i 2 . We see, first, that the perturbative 
contribution is significant and sometimes dominant and, second, that the modifications of DDIS 
factorisation are important, that is, the inhomogeneous term in the evolution equation plays a 
crucial role. 

It is informative to study the various contributions to dF^ /dh\Q 2 in some detail, since 
this is the quantity which mainly drives the behaviour of the diffractive gluon distribution. 
Recall that in the global analyses of inclusive DIS data, the small-x gluons are extracted using 
DGLAP evolution in which <9-F2/<9hiQ 2 ~ as g. As we have mentioned, in DDIS, the derivative 
dF^ 3 ^ /d\nQ 2 contains an additional contribution arising from the inhomogeneous term; see 
(fT2|) . In Fig. [7| we show the (3 dependence of the hiQ 2 derivative of F^^, for xjp = 0.003 and 
Q 2 = 15 GeV 2 , obtained in the recent MRW analysis [7] which takes [io — 1 GeV. From this 
figure, we see that the major contribution to dF^^ /d\nQ 2 comes from the first term in (|T2*j) . 
that is, the perturbative DGLAP contribution for fi > fiQ given by the dotted curve. The non- 
perturbative contribution coming from fi < fi is rather small for (3 < 0.5. The contribution 
corresponding to the differentiation of the upper limit of the /i 2 integral, that is, the second term 
in ()12j). is shown by the short dashed line. Clearly it is not negligible. The twist-four (F^ 3 ^) 
contribution is important, as expected, for large (3, and the contribution from the secondary 
Reggeon at this small value of xjp is practically invisible. In summary, we see that the slope 
for (3 < 0.2 is dominated by the perturbative DGLAP contribution, and for (3 > 0.8 by the 
twist-four F^ 3 ' contribution, while for (3 ~ 0.5 the inhomogeneous contribution is largest. In 
any complete analysis of DDIS data, it is clear that the inhomogeneous term must be taken 
into account. 

7 Discussion 

The QCD analysis of diffractive structure functions is subtle. To clarify the situation it may 
be helpful to pursue the analogy with the parton distributions of the photon. The evolution 
can be expressed in two alternative forms. In the first form, the DGLAP evolution (jSj) of the 
parton distributions of the photon includes an inhomogeneous term arising from the direct 
splitting of a point-like photon into a qq pair. This splitting may take place at any current 
scale. An alternative way to account for this effect is to treat the photon as a parton, and to 
consider homogeneous DGLAP equations which embrace the evolution of the gluon and quark 
distributions, together with that of the photon itself. For such a photon density we have the 
trivial distribution 7 = 8(1 — x), but now the inhomogeneous term is embodied in the enlarged 
set of homogeneous DGLAP equations for q, g together with 7. The advantage of the latter 
approach is that, in such a form, we have enlarged the subspace where the factorisation theorem 
applies. We can now factor off the matrix element of the hard subprocess i + j — > X where 
hj — <2S 9i or 7- Thus we have the possibility of the photon directly participating in the hard 
interaction. 

In the case of the diffractive parton distributions, we have a very similar situation to that 
of the parton densities of the photon. The DGLAP evolution contains an inhomogeneous 
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Figure 6: The perturbative (IP = G,S, and the GS interference) and non-perturbative contri- 
butions to F^ 3 \ for two sets of xjp, (3, and Q 2 values, found in the analysis of HERA DDIS 
data in Ref. [7]. The plots show the /i 2 dependence of the perturbative contributions; their 
integral over /i 2 is shown by the numbers in parentheses in the legend. 



14 



Xjp = 0.003, 



= 15 GeV 



Pert. dF^/dln Q 2 term 

Inhomogeneous term 

Non-perturbative 
Twist-four 

Secondary Reggeon 




Figure 7: The breakdown of the contributions to the description of the slope dF^ 3 ^ /d\nQ 2 as 
a function of (3, for Q 2 = 15 GeV 2 and xjp = 0.003 obtained in the MRW analysis [7] of the 
combined ZEUS and HI DDIS data, which uses MRST2001 NLO partons [16] to calculate the 
perturbative Pomeron flux. 
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contribution arising from the Pomeron-to-parton splitting, as given by the last term in (|10j). In 
the Q 2 — > oo limit, and at fixed xjp, due to the 1/Q 2 factor explicit in fjp, this may be treated as a 
power correction and safely neglected. Then the DDIS factorisation theorem is valid. However, 
for small xjp, the parton densities increase rapidly with increasing Q 2 , partly neutralising the 
power suppression of the inhomogeneous Pomeron-induced term. As a consequence we must 
retain the inhomogeneous term, and we are forced to reject the attractive postulate that DDIS 
factorisation [4] holds down to low scales. 10 

Since, in practice, it is easier to work with homogeneous DGLAP equations, we can proceed 
in one of two ways. First, just as in the photon case, we may treat the Pomeron as an extra 
parton, and enlarge the set of DGLAP equations to include fjp, as well as q° and g B . We denote 
the Pomeron density as fjp to be consistent with our previous notation. The extra terms P q jpfjp 
and Pgjpjjp in the homogeneous DGLAP equations for q D and g B involve the Pomeron-to-parton 
splitting functions. Moreover, just as in the photon case, to calculate the cross section we have 
to include the direct Pomeron coupling to the hard subprocess (see Section l4~2l for an example). 
In this approach the effective Pomeron distribution fjp is now described by the Pomeron flux 
fjp = fjp = c of (J7J), together with the corresponding flux fjp=s for the sea-quark component of 
the Pomeron, and the IP = GS interference term. Recall that the inclusive parton densities 
of the proton, xjpg and xjpS, which occur in these fluxes, satisfy their own DGLAP evolution 
equations. 

Equally well, the inhomogeneous DGLAP evolution of diffractive parton densities could be 
accomplished in an alternative way. We may sum up the parton densities given by a series of 
homogeneous DGLAP equations with different values of the starting scale \i. Suppressing the 
xjp and (3 dependence, the general structure is 

where aF^p;^) results from DGLAP evolution up to fi F starting from the input (fi; fi) = 
Pojp(P). The non-perturbative term, a D on _ pert , is exactly analogous to the procedure used in 
the Regge factorisation analysis described in Section [TJ but with «zp(0) fixed at the 'soft' value 
of 1.08 [6]. The important new ingredient is the pQCD contribution given by the second term 
on the right-hand side of (|17j) , which has contributions from the gluonic and sea-quark Pomeron 
(IP = G,S), and their interference {IP = GS). Setting ixp = Q, it is clear that (fTTj) satisfies 
the inhomogeneous evolution equation (fTOJ). 

In summary, we have shown how to obtain universal diffractive parton densities aP which 
can be used in the description of different diffractive processes. Of course, for diffractive 
production in high energy hadron-hadron collisions, we have to take care that the rapidity 
gaps are not populated by secondaries produced during the soft interaction of spectators. This 
well-known rapidity gap survival factor, usually denoted by S 2 [20], may be calculated from 
phenomenological models tuned to describe elastic and related 'soft' hadron-hadron processes 
[21,22], 

10 This postulate is the basis of the description of DDIS that is presented in Ref. [19]. We thank Paul Hoyer 
for clarifying discussions concerning this work. 
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In addition to the hard subprocesses originating from the collision of quarks and gluons from 
the resolved Pomeron, we must also include the contribution where the perturbative Pomeron 
directly participates in the hard interaction. In particular, the calculation in the Appendix 
of the cross section for diffractive open charm production may be treated as the coefficient 
function for 7*ZP — ► cc in the fixed flavour number scheme. 
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A Appendix: The gg and qq Pomeron in pQCD 

In this Appendix we compute the lowest-order Feynman diagrams which give the f3 dependence 
of the quark singlet and gluon starting distributions of the perturbative Pomeron, E ff '(/3, /i 2 ; fi 2 ) 
and (^(/^/z 2 ; fi 2 ). Here, the quark singlet distribution = vF + dF + s F + v! p + d lp + s IP , with 
all six quark densities of the Pomeron equal to each other. In Section lA.ll we compute Yf =G and 
(7 7P=G assuming that the QCD Pomeron is made from two t-channel gluons. Wherever possible, 
we check our results with those of Wusthoff [15] in the limit of massless quarks and with Ref. [23] 
in the case of massive quarks; see also the original work of Ref. [24] . Then in Section IA.2I we 
extend the formalism to calculate TF =S and cf p= assuming that the Pomeron is represented by 
a t-channel sea-quark-antiquark pair. Finally, in Section IA.3I we calculate YF =GS and ^ P=G5 
resulting from the interference between the gluonic Pomeron and the sea-quark Pomeron. As 
explained in the main text, the starting distributions, S iP (/5, /i 2 ; fi 2 ) and g^iP, fi 2 ; fi 2 ), can be 
thought of as Pomeron-to-parton splitting functions, P q jp and P g jp respectively. 

We will work in the leading logarithmic approximation (LLA) to derive the factorised form 
of (dU): 

<St.(^,Q 2 ) = / 9 /ip(W) F?((3,Q 2 ^ 2 ). (A.l) 

Higher order corrections (for example, to NLO accuracy) may be incorporated by considering 
proton and Pomeron parton densities satisfying NLO DGLAP evolution, by using the NLO 
coefficient functions to calculate F^, and by computing 0(as) corrections to the Pomeron-to- 
parton splitting functions. However, note that the factorised structure of (|A.1J) persists even 
at NLO (and beyond). 
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Figure 8: (a) Quark dipole and (b) effective gluon dipole interacting with the proton via a 
perturbative Pomeron represented by two t-channel gluons. 

A.l Two-gluon exchange (IP = G) 

First we consider the kinematics of the quark dipole shown in Fig. IHfa). We use a Sudakov 
decomposition of the momentum k of the off-shell quark, 



k = a q' + Pk V + k j 



(A.2) 



with q' = q + x B p, q' = = p , k]_ = —k\. The two outgoing components of the dipole have 
momenta 



(k 2 + m 2 f )/Q 2 
q — k = (1 — a) q + x B — p — kj 



1 — a 



(k? + m})/Q 2 
k = k + xjpp = aq + x B p + fej 



a 



(A.3) 
(A.4) 



where the on-shell conditions, (q — k) 2 = m 2 r and k 2 = mi, determine 



(k 2 + m 2 )/Q^ 

1 + J - 

1 — a 



13k = -X B | I ; ) • •'•;/' -''a 

The invariant mass of the qq system is given by 

M\ = (q + xpp) 



(k 2 + m 2 )/Q 

1 H 7Z 



a(l — a) 



(A.5) 



■2 , ,2 $ + m j 

a(l — a) 



(A.6) 



The kinematic limit occurs when a = 1/2, giving a maximum value for k 2 of Mj^/4 — mj. Since 
/3 = Q 2 /(Q 2 + M|), (|A~6|l can be written as 



a(l-a)Q 2 = P 



kf + mj 



P^ 2 



(A.7) 
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where the last equivalence specifies our choice of factorisation scale /x. The off-shell quark with 
momentum k in Fig. Ufa) has virtuality given by 



k 2 - m 2 = ^— ~ V, (Ai 

J 1 — a 



since a C 1 in the approximation of strongly-ordered transverse momentum, /i 2 <C Q 2 , to 
which we are working. 11 

A. 1.1 Quark dipole with a transversely polarised photon 

The differential 7*p cross section corresponding to Fig. |Sf a) is given by a fc r factorisation for- 
mula. It can be written in terms of photon wave functions k t ), describing the fluctuation 
of the photon into a quark-ant iquark dipole, convoluted over a and k t with a dipole cross 
section a, describing the interaction of the dipole with the proton via two-gluon exchange. The 
dipole factorisation formula for Fig. IBfa) with a transversely polarised photon is 

2 



N c f\ fdk 2 ^ 2 1 ^ fd% nT 

17 U 17 f J, LI 17 



167T ... , _., 

t=0 u / 7,h,A'=±l 



d 2 ^ „ T7 drj 



(A.9) 



dt 

where the qq dipole wave functions, ^1 h i, and the operation D are specified below. This formula 
for diffractive DIS may be compared to the corresponding result for inclusive DIS, 12 



/ 7,A,ft'=±l " " "* 

(A.10) 

The light-cone wave functions for the quark-antiquark dipole with a transversely polarised 
photon are [15, 23] 

, T ,7 (n u\_ 5h >~ h ' [( 1 ~ 2a ) h ~ 7] e 7 • fc t + <W m f h 

fetj " k 2 + m) + ajl - a)Q 2 ' (A - U) 

where 7, h, and ft/ denote the helicities of the photon, quark, and antiquark, respectively. Here, 
k t = k\ + i k 2 = (fc 1 , k 2 ) and the circular polarisation vectors of the photon are e 7 = (1,7)/ v2- 
The denominator of these wave functions is the virtuality of the off-shell quark with momentum 
k: 

k 2 t +m) + a(l - a)Q 2 = (1 - p)fx 2 + /3/i 2 = /i 2 ~ \k 2 - m)\. (A.12) 

Note that the wave functions are symmetric under a. — > (1 — a) and k t — > — k t , corresponding 
to q <-> g, that is, ()A.3|) <-> ()A.4|) . so we only need to sum over flavours in (jA.9|) and not over 
quarks and antiquarks separately. 

The four different permutations of the couplings of the two t-channel gluons to the two com- 
ponents of the quark dipole, shown in Fig. EJa), are obtained by simply shifting the argument 
of the wave function: 

mja^kM =2tt( q, k t ) - tt(a, k t + l t ) - (a, k t - l t ). (A.13) 

Actually, from (|A.7fl . fi 2 <§C Q 2 implies either a < 1 or (1 - a) < 1, but it is conventional to take the 



11 

former 



12 Note the extra factor 2 in (|A.10|I compared to (7) of [15]. 
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Figure 9: The four different permutations of the couplings of the two t-channel gluons to the 
two components of (a) the quark dipole and (b) the effective gluon dipole. 

We choose a basis where k t = k t (1, 0) and l t = l t (cos 0, sin 0) and neglect the xjpp components 
of the momenta. We work in the approximation of strongly-ordered transverse momenta, l t <C 
/i <C Q, and expand D^> in the limit U — > 0, only keeping the leading term proportional to I 2 . 
After performing the azimuthal integral, we find 



2tt 



27T 



D^l h ,(a,k t ,l t ) 



I 2 



x 



[k 2 + m 2 + a(l - a)Q 2 ] 6 
{4[a(l - a)Q 2 + m)\ 6 h> ^ [(1 - 2a)h - 7 ] e 7 • k t 

+2 SfM m f h [a(l - a)Q 2 + m) - k 2 ] } ; (A.14) 

cf. (21) of [15]. After changing variables from a and k 2 to (3 and fi 2 using (jA.7|) . we obtain 



7,h,/i'=±l 



2tt 



+mj m J+^3-0/i 2 J. (A.15) 



The other necessary part of the calculation is the cross section for qp —>■ qp, which may be 
obtained from that for the process qq — > gg with t-channel gluon exchange, shown in Fig. HUT a). 
We assume that k 2 ~ and neglect the p components of the momenta in these calculations. 
The qq — *> qq differential cross section is 



da \M\ 2 

W\ iqq ^ qq) = ^ 



(A.16) 



20 



(b) 




k 



fi,A u,B .p,B a, A 



k 




Figure 10: Cut diagrams giving the dipole cross sections for the two-gluon Pomeron: (a) 
qq — > qq and (b) gq — > gq. 

where s = 2 k -p = a Q 2 /xb — l^ 2 /xjp and t = —I 2 . The squared matrix element for qq — > qq is 



K9 4 



where the colour factor is 



C(qq -> qq) 



1 „ r.A .B-irrx. r .A + Bl _ ^ 



N 2 1 11 



9' 



(A.17) 



(A.18) 



In the high-energy limit, where only the terms to leading 0(\i\/s) are retained, that is, in the 
It — > limit, we have 

-^{qq -> qq) = -j—a s (l 2 )a s (fi 2 ), (A.19) 



dl? 



n 9 



where appropriate scales for as have been chosen corresponding to the two different vertices 
in Fig. [TOT a). The lower vertex in Fig. HUT a) may be considered as the first step of DGLAP 
evolution, which generates the unintegrated gluon distribution of the proton, f g (xjp,l 2 ,/i 2 ). 
Therefore, we obtain the cross section for qp — > qp by making the replacement 



2tt 



XlpP gq (X]p) 



xjp<1 



^2C F = ^-a s (l 2 )^f g (x IP ,l 2 ^ 2 ). 



(A.20) 



This replacement accounts for more complicated diagrams than Fig. ITUT a) which include the 
complete DGLAP evolution, leading to 



^(qp -> qp) = ^^-asi^fgixipJt,^ 



(A.21) 



Combining (jA.15|) and (jA.21|) we obtain 



j,h,h'=±l 



d 2 L 



7T 



da 

2 
t 



64 

,12 



/' J 

256vr 4 
9/i 12 



{...} 



2tt 2 dl 2 

— "s(m 2 ) -^-fgixpj 2 ,^ 2 ) 



[ as{^ 2 ) %ipg{xip, yu 2 ) ] 2 {• • •} , 



(A.22) 
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where g(xip,fi 2 ) is the integrated gluon distribution of the proton and {. . .} denotes everything 
inside the curly brackets in (jA.15|) . Strictly speaking, this last expression should be written in 
terms of the off-diagonal (or skewed) gluon distribution of the proton, since the left and right 
t-channel gluons in Fig. IHfa) carry different fractions of the proton momentum. At small Xp, 
and assuming that xjp g[xjp, /i 2 ) oc x]p X , then the off-diagonal gluon distribution is given by the 
diagonal distribution multiplied by an overall constant factor [12], 13 

2 2A+3 r(A + 5/2) 



R g (X) 



7T 



r(A + 4) 



flQ 2 



Changing the integration variables in ()A.9|h 

da dkf = d(3 du 2 - 

then integrating (|A.9J) over t, assuming a t dependence of the form exp(Bnt), gives 



(A.23) 



(A.24) 



TV 



R 9 N C 



T ' qq B D 167T 



dp 



d/i 2 /2 2 /Q 2 

2tt ^/I^I^JQ- 



f 



y,h,h'=±l 



d 2 L 



71 



The relation between the diffractive structure function and the 7*p cross section is 

F D( 3 ) = Q 2 g da^ 
kK 2 a era X]p d/3 

Thus, on inserting (|A.22|) into (jA.25|) we finally obtain 

Wd/J 2 1 1 T r> a s{H 2 ) 



7 D(3) 
T,qq 



E' 



/i 2 y/l - 4(3fi 2 /Q 2 XpBjo [ fi 



Rn 



xjpg{xjp,[i 2 



3/x 6 1 1 



(A.25) 
(A.26) 

(A.27) 



where again {. . .} denotes everything inside the curly brackets in (|A.15[) . Replacing X]/ e / by 
e 2 , this equation corresponds to diffractive open charm production; cf. (45) of [23]. 

We now take the limit of massless quarks, to/ — > where f = u,d, s, in the approximation 
where fi 2 <C Q 2 . In this approximation, we can replace the upper limit of the fi 2 integral, 
Q 2 /(4/3), by Q 2 . We also need to replace the lower limit by an infrared cutoff /i 2 ,. Then (jA.27|) 
becomes 



7 D(3) 
T,gg 



^d/x 2 1 
M 2 xjpB D _ 



R, 



o^sif 2 



xjpg{xjp,[i 2 



(e 2 )/3 3 (l-/?), 



(A.28) 



where (e 2 ) = G^/ e /)/ n / (with n f — 3), which coincides with (22) of [15]. We can write this 



expression as 



? D(3) 
T,gq 



d/i 2 



/jP =G (W)Ff= G (/3,/i^) 



where the 'Pomeron flux factor' is 14 

fjp=a{xjp\ ii 2 



xjpB 



D 



R, 



astjJ 2 



x IP g{x I p,iJ, 



(A.29) 



(A.30) 



13 This factor is not seen in calculating Fig. HUf a') since it is absent in the limit xjp — > [25]. 
14 Note that this definition differs slightly from that given in [7,8]. Since only the combination (|A.29|) matters, 
we are free to redistribute factors of fx and constants as we please. 
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and the 'Pomeron structure function' at a scale fi originating from a component of the Pomeron 
of size is 



F^= G (P, /i 2 ; /i 2 ) = (ej) PE F = G (P, /i 2 ; ^ 2 ) = (ej) /3 3 (1 - /?). (A.31) 

Recall that the notation IP = G is used to indicate that the perturbative Pomeron is represented 
by two i-channel gluons. The quark singlet density (|A.31|) can be taken as the initial condition 
at a scale /i 2 , together with the corresponding gluon density (|A.58J) . for DGLAP evolution up to 
Q 2 , giving JF= g ((3,Q 2 ;/j 2 ) and f =G {^, Q 2 ; /i 2 ); see Fig. El In this context, (3TF= G ((3^ 2 ; fi 2 ) 
can be regarded as the Pomeron-to-quark splitting function, P q jp=G{P)', see (JHJ). 



A. 1.2 Quark dipole with a longitudinally polarised photon 

The dipole factorisation formula for Fig. IHfa) with a longitudinally polarised photon is 

N c 



16tt 



t=o 



da f^rI2 e2 f a ™ 2 



h,/i'=±i 



7T dZ? 



(A.32) 



where the light-cone wave functions for the quark- ant iquark dipole with a longitudinally po- 
larised photon are [15,23] 



2a{l - a)Q 



8h,-h'- 



(A.33) 



k 2 + mj + a(l - a)Q 2 

Taking the limit If — > of the combination of wave functions (|A.13|) and performing the 
azimuthal integral gives 15 



Jo 271 



a(l — a)Q 2 + m 2 — k t 



[k 2 + m 2 + a(l - a)Q 2 ] 
Squaring this expression and summing over the quark helicities gives 



4a(l - a)Q 8 h) - h ' 



(A.34) 
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h,h'=±l 



2tt 



2tt 



4 32^ 
1 /i 6 Q 2 



P 2 [2(3 - 1 + 



2m 



2\ 2 



(A.35) 



The dipole cross section, da/dl 2 (qp — > qp), is the same as (jA.21|) . giving 



,D(3) 

~ Z^ C / /™ 2 , ,,2 

/ 



4/3 d/i 2 



/x 2 1 



£L /x 2 ^1 - 4/^ 2 /Q 2 Q 2 



x IP g{x I p,fi 



2 A3 



tf 3 / 2m 
^ 2/3 - 1 + 



(A.36) 

which corresponds to diffractive open charm production, cf. (44) of [23], after replacing Ylf e "f 
by e 2 . As before we take the limit of massless quarks in the approximation where fi 2 -C Q 2 , 
then 



,D(3) 



Q 2 (3 dal 



L ^ ^ 2 a era xjp d(3 



d/i 2 fi 2 



/i 2 Q 2 



(A.37) 



15 Note the extra factor 2 in (|A.34(1 compared to (21) of [15]. 
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where Ff =G ((3) = (ej) /? 3 (2/3 — l) 2 ; this result is a factor 2 different from (23) of [15], but 
is in agreement with [23]. Note that this contribution to F^ 3 ^ is twist-four due to the extra 
factor fi 2 /Q 2 with respect to (|A.29j) . Here we include only this main extra 1/Q 2 dependence, 
and omit the 0{ots) anomalous dimension arising from twist-four DGLAP evolution. 



A. 1.3 Gluon dipole with a transversely polarised photon 

Now consider the kinematics of the qqg system shown in Fig. Hfb). Although this diagram 
has an extra factor as with respect to the qq system shown in Fig. IHfa), it is known to be 
dominant at large Mx (small 0) due to t-channel spin-1 gluon exchange. Using a Sudakov 
parameterisation (jA.2|) of the momentum k of the off-shell gluon gives the momenta of the 
outgoing qq pair and gluon as 

™ p i J V T2_ 

q-k=(l-a)q'+ f 2 ^Y^f V ~ ^' (A ' 38) 

k 2 IQ 2 

k = k + xpp = aq' + x B — p + k ±J (A. 39) 

a 

respectively, where M qg - is the invariant mass of the qq system. Again the on-shell conditions, 
(q — k) 2 = M 2 q and k 2 = 0, determine 

/ k 2 + M 2 \ ( k 2 + aM 2 \ 

The invariant mass of the qqg system is given by 

M 2 x = (q + x IP p) 2 = ^- f. (A.41) 

a(l — a) 

The maximum value of k 2 occurs when a = (1 — M 2 ^/M x )/2, giving an upper limit for k 2 of 

Defining the (light-cone) fraction of the Pomeron's momentum carried by the gluon with mo- 
mentum k in Fig. Efb) to be 

k 2 /Q 2 * 



P= \x B -^-\ /xjp, (A.43) 

then the off-shell gluon with momentum k in Fig. IHJb) has virtuality 

k 2 = -k 2 /p = -ii 2 . (A.44) 

Referring to Fig. El the same kinematics hold if we now replace M q q by the invariant mass of 
all the emitted partons above the one labelled k in Fig. 01 

The qqg calculation is greatly simplified in the approximation M q q <C Q, in addition to 
assuming strongly-ordered transverse momenta, l t k t Q. In this limit, the kinematics of 
the (qq)g system is identical to the previously considered qq system of Fig. Efa) in the case 
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of massless quarks. The emitted (qq) pair is localised in impact parameter space, and forms 
an effective 'gluon' conjugate in colour to the emitted gluon. The (qq)g system can thus be 
considered as forming an effective gg dipole. This argument generalises if the (qq) pair is 
replaced by the emitted partons above the one labelled k in Fig. El all of which are strongly 
ordered in transverse momentum (see [9,11]). 

The dipole factorisation formula for Fig. IHfb) with a transversely polarised photon is then 



dt 



N 2 c 



16tt 



t=o 



da 



dk 2 t ^ 2 
f 



m,n=l,2 



d 2 Z 



7T 



dl 2 



(A.45) 



Again there are four different permutations of the couplings of the two gluons to the two 
components of the effective gluon dipole, shown in Fig. HU^b), which are obtained by shifting 
the argument of the wave function as in (jA~T3jl . The light -cone wave functions for the effective 
gluon dipole with a transversely polarised photon are [15] 



<H mn (a,k t ) 



k 2 5 mn - 2k™k? 



m, n 



1,2. 



y/a(l -a)Q 2 k 2 t + a(l - a)Q 2 
Taking the limit If — > of D^f mn and performing the azimuthal integral gives 



j^D^ mn (a,k t ,l t ) = l 



2k 2 t 



3ft(l - a)Q 2 + k 2 



y/a(l-a)Q 2 [k 2 + ft(l - a)Q 2 f V " k 2 
cf. (24) of [15]. Squaring this expression and summing over the indices m, n 



(A.46) 



(A.47) 
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The dipole cross section for gp — > gp is obtained from the scattering process gq 
i-channel gluon exchange, shown in Fig. fTUT b) . Here, the squared matrix element is 



1, 2 gives 

(A.48) 
gq with 



\M\ 2 = ~C^ Tr[^ 7 (^ - hh 5 } d^(k,p) d vp (k + l ± ,p) (-g aj ) (-g ps ) 



(A.49) 



x [(2k + l±) a g» u - (k + 2l i _Yg ua + (l ± - k) u g pa ] 

x [(2k + l ± y g pa - (k + 2l 1 _Yg pf3 + (l ± - k) p g uf3 ] 

where the transverse polarisations of the incoming and outgoing gluons are summed in a light- 
cone gauge, 



d^(k,p) = -g^a + 



kfiPa + Pfj,k a 



and where the colour factor is 
C(gq -> gq) 



N C (N ( 



p ■ k 

f ABC f ABD Ti[t c t D ] 



1 

2' 



(A.50) 



(A.51) 



In the high-energy limit, where only the terms to leading 0(s/\t\) are retained, that is, in the 
l t — > limit, 

■^(gq -> gq) = ji^ s (i 2 )a s (fi 2 ) 



(A.52) 



n 
25 



where appropriate scales for a$ have been chosen corresponding to the two different vertices 
in Fig. ITUT b). As before, we obtain the cross section for gp —* gp by making the replacement 
(|A.20p . which gives 

^(gp -> gp) = ji^Y<x s (ii 2 )f g (xip, ll (A.53) 

Inserting (|A.48|) and (|A.53|) into (|A.45j) . and accounting for the skewed effect, the change of 
integration variables from (a, fc 2 ) to (/3,/x 2 ), and the t dependence, we obtain 

da^ p 4vr 2 a cm f Q2 dp, 2 



dp g 2 



FKwl^ 2 ^^**^] 2 E e /I( 1 -/ 3 ) 2 (l + 2/3) 2 i (A.54) 



We must account for the fact that the off-shell gluon with momentum k in Fig. HJb) does 
not interact directly with the photon, but first splits into a quark-ant iquark pair forming the 
effective 'gluon' of the dipole. To do this, we replace j3 — > /3' in the previous formula and 
include the DGLAP splitting for g — > qq, that is, 

A °T,(,4)g _ a s (Q 2 ) ( Q 2 \ f'dff' D (f}\ d<4> 



w 2, ln w ]„ f p " \w) (A - 55) 

Putting everything together, we finally obtain 

F D(3) Q 2 P d(T T,( q g)g _ f^dfl 2 2 jp =G , 2 , 

where 

jf= G (/3,Q 2 ;^ 2 ) = <e 2 }/?E p = G (/3,QV 2 ) 

= 2£4^<« (f ) (f ) A (A.57) 

with 

/3^= G (/3', /i 2 ; /i 2 ) = A (1 - /5') 2 (1 + 2/?) 2 . (A.58) 

lb 

Note the extra factor 2 compared to (25) of [15] (this was corrected in a later paper [26]). We 
take (|A.58J) at a scale /i 2 , which can be interpreted as the Pomeron-to-gluon splitting function 
P g jP=G ©, as the initial condition for DGLAP evolution up to Q 2 . 



A. 2 Two-quark exchange (IP = S) 

We now calculate the lowest-order Feynman diagrams for the case in which Pomeron exchange is 
represented by a t-channel sea-quark-antiquark pair, rather than two gluons. The quark dipole 
and effective gluon dipole interacting with this two-quark Pomeron are shown in Fig. ^2 The 
light-cone wave functions of the photon, ^f(a, k t ), are the same as those given in Section fA. II 
The two different permutations of the couplings of the two sea quarks to the two components 
of the dipole, shown in Fig. axe obtained with 

D*(a,fe t ) = 2*(a,fc t ), (A.59) 
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p' P 

Figure 11: (a) Quark dipole and (b) effective gluon dipole interacting with the proton via a 
perturbative Pomeron represented by two t-channel sea quarks. 
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Figure 12: The two different permutations of the couplings of the two i-channel sea quarks to 
the two components of (a) the quark dipole and (b) the effective gluon dipole. 
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Figure 13: Cut diagrams giving the dipole cross sections for the two-quark Pomeron: (a) 
qq -> gg and (b) gq -> qg. 

that is, there are no terms with a shifted argument, unlike for the two-gluon Pomeron. Since 
there is no It dependence here, the integrals over l t in the dipole factorisation formulae can be 
done immediately. 

A. 2.1 Quark dipole with a transversely polarised photon 

The dipole factorisation formula for Fig. HIT a) with a transversely polarised photon is 



dt 



N, 



c l\_ fdh? 



da 

lQn In J 2tt 



\ E w^. 



(A.60) 



fe,h',7=±l 



where 



E i^(«>M 2 = ^{[(i-/V-mj] ( 1 - 2 / 3 ^)+ m /|- ( A - 61 ) 

7,/i,h'=±l ^ ^ V V / J 

The dipole cross section for qp — > is obtained from the scattering process gg — > with 
i-channel sea-quark exchange, shown in Fig. 113( a) . Here, the squared matrix element is 



\M\ 2 = \C^ Tr[ 7 ^ 7 7±7WW M* + 4*,(p - 'JL. 



where the colour factor is 



1 1 

C(gg-gg) = ^Tr[tWt B ] = ^. 



(A.62) 



(A.63) 



In the high-energy limit, where only the terms to leading 0(s/\i\) are retained, 



da 
df 



(qq -> gg) = 4^^|a5(^)«s(^ 2 )- 



(A.64) 



Note that this expression is suppressed by an extra factor 0(\i\/s) compared to the t-channel 
gluon exchange processes considered in Section |A~TI where s ~ /i 2 /xjp and t = —I 2 . We obtain 
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the cross section for qp — > #p by making the replacement 



2tt 



XlpPqq(X]p) 



271 07T 



(A.65) 



where f q (xip,l 2 , fi 2 ) is the unintegrated quark distribution of the proton. This replacement 
gives 



a(qp -> 



167T 2 as{fi 2 
~~9 ^ 



/■ M d/ 2 

/ -j£-fq(X]P,%,tl 2 ) 

Jo n 



16n 2 as{fi 2 
~9 



-xjpq{xjp, /j 2 ), 



(A.66) 



where q(xjp,fi 2 ) is the integrated quark distribution of the proton. Again, we should use the 
skewed quark distributions, which gives rise to an extra factor [12] 



R q (X) 



2 2A+3 r(A + 5/2) 
V5F T(A + 3) ' 



(A.67) 



assuming that xjp q(xjp, /i 2 ) oc x F x at small xjp. The final formula for massive quarks is 



7 D(3) 
T,qq 



1 



1 



4/3 dyU 2 

4^ s /l^lf3 l i 2 IQ 2 x 1P B D 



xpqlxp, fj 2 



16(3 
2V 



{...} 



(A.68) 

where {. . .} denotes everything inside the curly brackets in (|A.61J) . Replacing 2^/ e / ^y e 2 
and q(xip,^ 2 ) by c(x]p,fi 2 ) in ()A.68J) gives a formula for diffractive open charm production. 
Alternatively, assuming light quark flavour symmetry, q(xjp, fi 2 ) = S(xjp, fi 2 )/(2rif) with rij = 3, 
and making the approximation fi 2 -C Q 2 we obtain 



P D(3) 
r T,qq 



dfi 2 
'1% » 

where the 'Pomeron flux factor' is 16 

fip=s(xip; fi 2 ) 



, fiP = s(xiP^ 2 )F^= s ((3^ 2 ^ 2 ), 



Rg <*s(p ) X]p s(x I p,fx 2 ) 
fi 



(A.69) 



XipP 



D 



(A.70) 



and the 'Pomeron structure function' at a scale fi originating from a component of the Pomeron 
of size 1/ fi is 



(/3,^/i 2 ) = <^> / 9E M C9 )A i a ;A* a ) = <ej>^ /?(!-/?) 



(A.71) 



Recall that the notation JP = 5 is used to indicate that the perturbative Pomeron is represented 
by two t-channel sea quarks. Again, (|A.71j) can be taken as the initial condition for DGLAP 
evolution and it can be regarded as the Pomeron-to-quark splitting function, P q jp=s', see (JT3J). 

16 Just as in l|A.30[) . we use a slightly different definition from that given in [7,8]. Since only the combination 
(|A.69(1 matters, we are free to redistribute factors of \i and constants as we please. 
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A. 2. 2 Quark dipole with a longitudinally polarised photon 

The dipole factorisation formula for Fig. ITTT a) with a longitudinally polarised photon is 

N c 

■L an I— . 



dt 



t=0 



16n 



/W^E^ E 

J ° J f h,h'=±l 



0|2 .2 



where 



h,/i'=±l 

leading, in the case of massive quarks, to 



v 2 Q 



r D(3) 
L,qq 



4/3 d/i 2 



1 



^ 1 



/i 2 y/i - A(3^/Q 2 Q 2 xjpB D 



R, 



a 5 (/i 2 ) 



l 2 



x ]P q{x I p,iJ J 



and in the limit of massless quarks, in the approximation |i 2 C Q 2 , to 



,D(3) 



^a^xjp d/3 



JP=S/ 



Q 2 



where Ff =s (f3) = (16/81) (e 2 ) /3 3 . Note that this contribution to F 2 D( - 3 ^ is twist-four due to the 
extra factor fi 2 /Q 2 with respect to (|A.69|) . and hence we do not perform leading-twist DGLAP 
evolution. 



(A.72) 



(A.73) 



|/? 3 , (A.74) 



(A.75) 



A. 2.3 Gluon dipole with a transversely polarised photon 

The dipole factorisation formula for Fig. [TTT b) with a transversely polarised photon is 

H_ 
2vr 



dt 



1 da 1 1 



t=o 



167T 



E2 
6/ O^em 



E l«* 



mn 1 2 « 2 
0" , 



m, n=l, 2 



where 



m,n=l,2 

The squared matrix element for gq — > gg, shown in Fig. IT3T b). is 



|.M| 2 = ^ Tr[ 7 7±7^77±7^(^ + W] MM rf P> ~ J±. *). 



where the colour factor is 



J TrftVt B t B l = - 



The dipole cross section is then 



^(gq -> gy) - , 2 



1 32vr XiP 2 2 



-a s {l 2 t )a s {n 2 ), 



(A.76) 



(A.77) 



(A.78) 



(A.79) 



(A.80) 
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so, making the replacement (|A.65|) . the gp — > qp cross section is 

u \ 27r 2 a 5 (/i 2 ) r' 2 d/ 2 ^ 2 2 2ir 2 a s (n 2 ) , 2 , . 

a(gp^qp) = — — / > f q {xjp,l t ,n ) = — —xpSixp, /i ), (A.81) 

6 ^ Jo n o fi 

assuming light quark flavour symmetry, q(x]p,fi 2 ) = S(xp, fi 2 )/(2nf) with n/ = 3. Again, we 
account for the g — > qq splitting using (|A.55j) . Putting everything together, we finally obtain 

J til V 

where 

if = s (/3, Q 2 | pi') = (eJ)fiTr- s {13,qi>;p?) 

= 2£4^rMf ^ S{p '' M ' (A - 83) 

with 

P'g JP = s (P'^ 2 ^ 2 ) = 1 -(l-P') 2 - (A-84) 

We take ()A.84|) at a scale /x 2 , which can be interpreted as the Pomeron-to-gluon splitting 
function P g jp=s flEl), as the initial condition for DGLAP evolution up to Q 2 . 



A. 3 Interference between two-gluon and two-quark exchange (IP = 

GS) 

We must account for interference between two-gluon and two-quark exchange, that is, interfer- 
ence between Fig. [HJ^a) and Fig. ITTlf a). and between Fig. IBfb) and Fig. ITTlf b). We label these 
contributions using the notation IP = GS. 



A. 3.1 Quark dipole with a transversely polarised photon 



The dipole factorisation formula for interference between Fig. Eta) and Fig. [TIT a), with a 
transversely polarised photon, is 



dt 



N, 



c 



t=0 



16n 



da 



d 2 L 



/ y,h,h'=±l 

Proceeding as before, the final result for massive quarks is 



7T 



\IP=S 



TP=G 



(A.85) 
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Q 2 



. (A.86) 
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In the limit of massless quarks, in the approximation /i 2 <C Q 2 , this expression reduces to 

" 2 d/i 2 



7 D(3) 
T,qq 



where the 'Pomeron flux factor' is 



2 /iP =G 5(^;// 2 )i^ =G5 (/3,/i 2 ;/i 2 ), 

C 2\ 1 2 



(A.87) 



(A.88) 



and the 'Pomeron structure function' at a scale /x originating from a component of the Pomeron 
of size is 



(A.89) 



Again, (|A.89J) can be taken as the initial condition for DGLAP evolution and it can be regarded 
as the Pomeron-to-quark splitting function, P q jp=GS\ see (HSj)- 



A. 3. 2 Quark dipole with a longitudinally polarised photon 

The dipole factorisation formula for interference between Fig. Efa) and Fig. flTT a). with a 
longitudinally polarised photon, is 

J f h,h'=±l 



dt 



t=0 



dl 2 
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IP=G 



(A.90) 



For the case of massive quarks, this leads to 



j,D(3) 
L,qq 



/m? 
/ J — 



d/i 2 



/i 2 R g R q 



"± /i 2 y/l - 4f3fi 2 /Q 2 Q 2 xjpBjj 
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2 xpgi^xp, /i 2 ) xjpq{xjp, /i 2 
8/3 3 



2/3-1 + 



2m 2 



9 V ^ 

and in the limit of massless quarks, in the approximation |i 2 C Q 2 , we obtain 



? D(3) 
L,qq 



Q 2 (3 d< 



TP_ 
qq 



Air 2 a em xjp dp 



d/i 2 /i 2 
~W Q 



<IP=GS 



(A.91) 



(A.92) 



where Ff= GS ((3) = (4/9) (ej) /3 3 (2/3 - 1). Note that this contribution to F 2 D(3) is twist-four due 
to the extra factor /i 2 /Q 2 with respect to (|A.87j> . and hence we do not perform leading-twist 
DGLAP evolution. 



A. 3. 3 Gluon dipole with a transversely polarised photon 



The dipole factorisation formula for interference between Fig. E^b) and Fig. fTTT b). with a 
transversely polarised photon, is 

dlV 



dt 



N 2 -l f\ fdk 2 
( 1 da 1 



167T 



t=0 







^ e "f ' ! 1-111 
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2 * Y' 



£ 2 



TO,n=l,2 



IP=G 



(A.93) 
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Again, we account for the g — > qq splitting using (jA.55j) . Putting everything together, we 
finally obtain 

F ?, ( &) 9 = /f ^ f^osixjp; /. 2 ) Ff= GS (f3, Q 2 ; (A.94) 

where 

JT cs (/3,QV 2 ) = (e^E^OS.QV 2 ) 

= 2 E4^to (f ) ? 0n» (I) tr<r«*V,M, (a.95) 

with 

^/^(/^V) = i (1 -/?') 2 (l + 2/3'). (A.96) 

We take ()A.96|) at a scale /i 2 , which can be interpreted as the Pomeron-to-gluon splitting 
function P g jp =GS (JTHJ), as the initial condition for DGLAP evolution up to Q 2 . 
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